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Locally inertial coordinates are constructed by carrying Riemann normal coordinates on a codi- 
mension two spacelike surface along the geodesies normal to it. Since the normal tangents are 
labelled by components with respect to a null basis, these coordinates are referred to as null normal 
coordinates. They are convenient in the study of local causal horizons. As an application, the 
coordinate system is used to specify a vector field that satisfies the Killing equation approximately 
in a small region and the Killing identity exactly on a single null geodesic. We also construct a 
vector field on a surface, starting from a vector at a given point on the surface. This construction 
may be regarded as a generalisation of the notion of Fermi- Walker transport. 
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Different problems call for different coordinate systems. For example, a fluid whose flow lines are hypcrsurfacc- 
orthogonal geodesies is suitably described in terms of Gaussian coordinates. These are defined as follows: On a 
(portion of a) spatial hypersurface f2 some arbitrary coordinates are placed. A point r in the neighborhood of Q lies 
on a unique timclike geodesic intersecting SI orthogonally, at q say. Then the Gaussian coordinates of r consist of the 
proper time from q to r along the geodesic, together with the spatial coordinates of q on SI. 

Null normal coordinates (NNCs) are defined in a similar way, but starting from a codimension two spatial surface 
E. At each point on E, the tangent space orthogonal to E may be spanned by a pair of null vectors. A point r in the 
£h , neighborhood of E again lies on a unique geodesic intersecting E orthogonally, at q say. At this point q, the tangent 
^j' to the geodesic will be a linear combination of the pair of null vectors. If the geodesic is parametrized such that r lies 
at unit afiine parameter, then the NNCs of r are given by the coefficients of this linear combination, together with 
' arbitrary coordinates of q on E. Accordingly, every null geodesic orthogonal to E is a coordinate curve. 

NNCs were introduced in [l|, in the context of black hole horizons, where E was chosen to be a cross section 
of the horizon. However, from the definition above it is clear that NNCs may be constructed off any codimension 
two spatial surface in an arbitrary spacctimc, as long as the surface can be covered with a single coordinate chart 
, and the orthogonal geodesies do not intersect. Note also that the orthogonal null geodesies generate a pair of null 
hypersurfaces, each cut in half by the surface. This calls to mind the "lightsheets" of the covariant entropy conjecture 
put forward in Q. 

Furthermore, the null surfaces to the past of E form part of the boundary of its past. Therefore, if E is very small, 
\ the null surfaces are part of a local causal horizon. Such local horizons have been employed to address the question 
whether gravitational field equations are a manifestation of the Clausius relation of thermodynamics 0, @ • 

In scction|TT]we will construct NNCs starting from a particular kind of small surface E, namely one that is generated 
by geodesies. More precisely, E is the image under the exponential map of a small part of a spatial, codimension two 
' tangent space of a spacetime point p. Starting from such a surface has the added advantage that the NNCs can be 
made to be locally inertial, i.e. at p the coordinate vectors are orthonormal 1 and the Christoffel symbols vanish. It 
will also be shown that the coordinates we introduce are, in some sense, the simplest possible to have the properties 
of being null normal and locally inertial. 

The construction of the coordinate system requires a pair of null vector fields orthogonal to the surface E. While in 
section |TI] these are specified with the aid of a Riemann normal coordinate system, such vector fields may equivalently 
be defined in a coordinate independent manner, via a generalisation of the notion of Fermi- Walker transport. We 
introduce this notion in the appendix. Starting from a given pair of null vectors at a point and transporting according 
the prescription given in the appendix uniquely retrieves the vector fields constructed in section [TI] 

An important ingredient in the conjectured thermodynamic nature of gravity is a collection of local, approximate 
Killing vector fields. The components of these vector fields in locally inertial coordinates are such that boosts would 
be generated if the coordinates were globally inertial. Then the Killing equation is satisfied in the neighbourhood of 
a spacetime point, up to possible corrections linear in the locally inertial coordinates. 
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1 See the remark below Eqn. H14I I for a qualification of this statement. 
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However, this does not specify the vector fields uniquely. Using the coordinate system of section [TT1 we will show 
in section I Jill that it is possible to further specify the fields such that they satisfy the Killing identity 2 exactly when 
restricted to a single geodesic. This additional specification also has the consequence that the Killing equation in a 
neighbourhood of a spacetime point holds up to corrections quadratic in the locally inertial coordinates, and holds 
exactly on the geodesic. 

Since the study of local horizons of Q makes use of vector fields satisfying these properties, the calculations reported 
in the present paper may be regarded as providing more explicit detail to the proof that such fields exist, as compared 
with the more sketchy treatment given in NNCs arc particularly well adapted to that study since the generators 
of the local horizon simply become coordinate curves, whilst further specialising the coordinates to be locally inertial 
naturally allows for the construction of an approximate boost generator. It is hoped however that locally inertial 
NNCs and the associated approximate Killing vectors may be of some use in other contexts as well. 



II. CONSTRUCTION OF LOCALLY INERTIAL NULL NORMAL COORDINATES 



At a point p of an n-dimensional spacetime, a basis {l a , k a , e A } , A = 2, . . . ,n — 1, of the tangent space is introduced, 
where {e A } is an orthonormal set of spatial vectors, and l a and k a are null vectors orthogonal to e a A . Furthermore, 



the null vectors are scaled so that their inner product equals minus one . In summary, 



and 



k a k — l a l — k a e A — l>a?-A — 0, ^ak 
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(1) 



(2) 



We assign the Riemann normal coordinates (RNCs, see @) {y a } = {u 7 v, y A } to a point in the neighbourhood of p 
if it lies at unit affine parameter on the geodesic through p with tangent vector at p given by u l a + v k a + y A e A . The 
components of the inverse metric are expanded as 



9 al3 = rf t3 + \R a f u y^ + 0(y 3 



(3) 



Here, O (y 3 ) denotes terms proportional to at least third powers of RNCs and r] a P is the flat metric in double null 
coordinates, namely 



a/3 



A U A- 



(4) 



Furthermore, throughout the text Riemann components are understood to be evaluated at the origin only 4 . 

Let E be the (n — 2)-dimcnsional spatial surface defined by the equations u = v = 0, i.e. E is generated by geodesies 
with tangent at p of the form y A e a A . We now look for null vector fields spanning the planes orthogonal to E; they 
will serve as coordinate vectors of the NNC system. Note that the parallel transport of e.g. k a along a generator of 
S does not in general remain orthogonal to £ in a curved spacetime, while the orthogonal projection of the parallel 
transport of k a is not guaranteed to be a null vector. On the other hand, the gradients of u and v are orthogonal 
everywhere. However, they are null only at p, with —V a u(p) = k a and — V a v(p) = l a . Therefore, we will look for 
linear combinations of the gradients of u and v that are null everywhere on £. Solving a quadratic equation and 
choosing the appropriate sign shows that the vectors 



K a = -V a it + KV a v, L a = -\7 a v + LS7 a u 
are null everywhere on E, provided that K and L take the values 



K = (g™)" 1 
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2 For a true Killing vector, the Killing identity V a Vj,5c = R d abdd obtains by taking covariant derivatives of the Killing equation 

V (a 4) = 0. 

3 We use the sign conventions of [f|, i.e. the metric signature is !- + ••■ and the Riemann tensor is defined by 2V[ a Vj,]o; c = R^^lu^. 

Lower case Latin indices are used for abstract index notation, whilst Greek indices denote coordinate components. Finally, upper case 
Latin indices denote coordinate components on a codimension two spatial surface. 

4 with the exception of equation 11331 . 
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and 



{g uuyl 



gUV + ^ g Uvf_gUU gl 



(7) 



Rescaling 5 the vectors so as to make their inner product equal to minus one, we obtain that for all q £ E the vectors 

(8) 



k a (q) = (-L b K b ) 1/2 K\ l a (q) = (-L b K b ) 1/2 L a 



are null, orthogonal to E, satisfy l a (q)k a (q) = —1 and at p reduce to the coordinate vectors of u and v. Anticipating 
the introduction of the NNC system we denote the spatial RNCs on E by 



x A = y A (q). 

The RNC components of the null vectors §8§ are expanded in terms of the spatial coordinates ^ as 
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k V = 1 - -RuAvB X A X B + O (x 3 ) , l U = 1 - -RuAvB X A X B + O (x 3 ) 
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k u = -\r vAv b x A x B + O (x 3 ) , l v = -\r u auB x a x b + O (x 3 ) 
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(9) 



(10) 



(11) 



k° = -RcAvB X A X B + O (x 3 ) , l C 
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= TiRcAuB X A X B + O (x 



(12) 



To introduce NNCs {x a } = { U, V, x A } , we note that a point r in the neighbourhood of p lies on a unique geodesic 
intersecting E orthogonally, at q say. Let the tangent to that geodesic at q be given by Ul a + Vk a when r lies at unit 
affine parameter (, sec figure 1. The NNCs of r are then defined by 



or in other words 



x a (r) = U6Z + V6Z + y A (q)5% 



x°(r) = U, x 1 ^) = V, x A {r) = y A {q) = x J 



(13) 



(14) 



These comprise a locally inertial coordinate system at p, as we will now show (Strictly speaking, it is the coordinate 
system {(U + V)/\/2, (U — V)/\/2, x A } which is locally inertial in the usual sense.). 

To find the lowest order terms of the metric components in the NNCs (fTl)) , we establish the coordinate transfor- 
mation between RNCs and NNCs by solving the geodesic equation perturbatively. Namely, the RNCs y a (r) of the 
point r, lying at affine parameter ( = 1 on the geodesic intersecting E orthogonally at q, are the solutions y a (q, 1) to 
the equation 



with initial conditions 



and 



d 2 v a dv 13 dyT 

-^(q, C) + T%,{q, C)^fe 0-^(9, = 0, 

u (q, 0) = v {q, 0) - 0, y A (q, 0) = y A {q) - x A 

dy Q , 



( q ,0) = Ul a (q)+Vk a (q). 



(15) 



(16) 



(17) 



Using the expressions (fTTJI) — (JT3J) for k a (q), l a (q) and expanding the Christoffel symbols of the RNCs in terms of the 
Riemann tensor Q , we solve this equation to obtain the RNCs at r in terms of the NNCs at r and the RNC components 



Explicitly, the rescaling factor is given by (— Li,K b ) 



2 \(g uv ) 2 -g uu g vv \ g uv + J{g uv ? -g uu g vv /(g uu g vv ) 
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qr, 

Ul a (q) + Vk°(q)J 
l a (q)' 




k°(q) 



Figure 1: The geodesic surface E in 3-dimensional spacetime (imagined to be a line segment perpendicular to the page). The 
null vectors k a and l a are orthogonal to E everywhere. The point p £ E is at the origin of the coordinate system and V is the 
null geodesic through p with tangent k a . At q £ E the orthogonal geodesic qr has a tangent Ul a (q) + Vk a (q) such that r lies 
at unit affine parameter The null normal coordinates (NNCs) of r are given by {x a (r)} = {U,V,y A (q)}, where y A (q) are 
Riemann normal coordinates of q on E. 



of the Riemann tensor at p. However, the result shows that to first order both coordinate systems coincide, so that 
dy a / dx@ (p) = 5p. The tensor transformation law then tells us that at p the NNC components of the Riemann tensor 
coincide with those in RNCs. Thus, with 



u,v 



the coordinate transformation may be written as 

1 



v(r) = V + -RuABi x i i' + -Ruvbi Vx B x l + O (x 
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u(r) =U+ -RvABi x A x M x l + -RvuBi Ux"x l + O (x 4 ) , 
6 3 

y C(r) = x c _ }- RcAm x A x B x l - ^R CjBi x J x B x l + O (x 4 ) . 
Applying the tensor transformation law to ^ then obtains the NNC components of the metric, to wit 

guv - -1 + ^Ruvuv VU + O (x 3 ) , 



9uu = ^Ruvvu V 2 + (x 3 ) 



gvv = ^Rvuuv U'+O (x 3 ) , 



(18) 

(19) 
(20) 
(21) 

(22) 
(23) 
(24) 



9ab = Sab + Razjb x % x 3 + ^R A {%C)b x l x c + -Racdb x c x d + O (x 3 ) (25) 



9UA = T^RuviA Vx l + ^Ruvba Vx B + O (x 3 ) , (26) 
9va = ^RvmA Ux i + ^Rvuba Ux B + O (x 3 ) . (27) 



A. Maximally simple NNCs 

Manifestly, the coordinate system defined by (QJJ is locally inertial and null normal as defined in the introduction, 
but not uniquely so. An infinite number of coordinate systems with these properties exists, even if the coordinate 
basis {l a ,k a ,e A } is kept fixed at p. Since no natural scale is associated to null geodesies, we may rescale the null 
vectors of Eqn. (JHJ through 

~k a (q) = f(q) k a (q), l a (q) = jL- l a {q) (28) 

where the rescaling factor f(q) is a smooth function on S, chosen such that f(p) = 1. The rescaled vectors are still 
null, orthogonal to S and satisfy l a k a = —1. They may again be used to construct NNCs, by the same procedure as 
outlined above. However, if d A f(p) 7^ 0, these coordinates will no longer be locally inertial. Therefore, for a given 
basis {l a , k a , e A } at p, the class of locally inertial NNCs obtains from the one defined above by restricting to rescaling 
factors of the form 

f(q) = l + f AB x A x B + O (x 3 ) , Jab € R. (29) 

Without loss of generality, the coefficients Jab may be taken to be symmetric in AB. 

Such a rescaling has no effect on the lowest order terms of the metric components except in giA, for which one finds 

9UA = 9UA - Vab Vx b + O (x 3 ) , g VA = g VA + 2f AB Ux B + O (x 3 ) . (30) 

One may wish to simplify the form of guA, by choosing a rescaling factor with appropriate coefficients Jab- However, 
the coefficient of the term proportional to Vx B in gjj A is given by i Rjjvba- This is antisymmetric in AB, so the 
term proportional to Vx B cannot be removed by a rescaling of the form (|29p . Similar remarks of course hold for 
gvA- Accordingly, the locally inertial NNC system defined by (Q3J) is the one in which the expression for the metric 
components is as simple as possible. 



III. CONSTRUCTION OF A LOCAL, APPROXIMATE KILLING VECTOR 

In a neighbourhood of a spacetime point p, we define a vector field £ a by taking its components in the NNCs (fH)) 
to be a power series 

T = V&v -U5Z + (x 2 ) , (31) 

such that the O (a; 2 ) terms vanish in Minkowski space. In Minkowski space, this defines the Killing vector that 
generates boosts in the UV plane, with fixed point at p. In curved spacetime no Killing vectors are guaranteed to 
exist, but p is still a fixed point of £ a at which the Killing equation is satisfied, and furthermore V( Q ,^) = O (x). 

We will impose a restriction on £° that determines many of the higher order terms in its expansion. For a true 
Killing vector, combining different permutations of the covariant derivative of the Killing equation V( a £(,) = results 
in the Killing identity 

V r V s („ — R m rsa£,m = 0. (32) 

For approximate Killing vectors of the form (f3"Tj) . the LHS of (|32l) is at best linear in the coordinates when evaluated 
in the neighbourhood of a point in a general spacetime. However, £ a may be chosen so as to make the LHS of (|3"2")l 
vanish when restricted to a single curve, as we will now show. 
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We begin by writing out the LHS of (|32j) in terms of partial derivatives and Christoffel symbols, as 

V,V CT £a - R% aa ^ = d p d^ a + {-d p Y^ a + 2Y" p( X a) J) ^ - (rV^k + 2r \( P d ^) ~ & poc.tr ( 33 ) 

This is an equation valid in any coordinate system and it is the only equation in this paper where Riemann components 
and Christoffel symbols take values in a spacetime neighbourhood, not just at p. Next we write the covariant 
components of £ a in NNCs as the power series 

£ tt = U5l - V5 u a + ^C paa x p x' 7 + -D^sSx"^ + -E papva x<>x a x»x v + O (x 5 ) , (34) 

where the coefficients C paai D papa etc. are symmetric in all but their last index. Lastly, let T be the null geodesic 
through p with tangent k a . We will evaluate the NNC components of the Killing identity term ([33]) on T, using the 
approximate Killing vector (|3"4"|) . Since Y is simply a coordinate curve, with points on Y having NNCs (0, V, 0, 0, ...), 
the NNC components of (|33|) on T reduce to a power series in V. Making use of the fact that at p the coordinate 
system is locally inertial, we find 

(VpV CT £ Q - R M paat^lr = C p™ + {Dvpoo + (paa) - (apa) - (aap) - d p Y Vaa - Rvpaa) V + (V 2 ) . (35) 

As in the rest of the text (with the exception of ([22])) the Riemann components and Christoffel symbols on the RHS 
of (J2H) are evaluated at p. Furthermore, we defined r a ^ 7 = g a gY s ^ and for brevity we introduced the notation 

(paa) = 2d v T [V p X ] - (36) 

Now we are ready to impose that the Killing identity hold on V, which amounts to putting the RHS of (f3"B"|) to zero, 
order by order. At zeroth order this requires 

C paa = 0, (37) 

and as a consequence the quadratic term in £ a vanishes, i.e. = VSy — US^ + O (x 3 ) . A further consequence is that 
the Killing equation is satisfied in a neighbourhood of p up to terms quadratic in the coordinates, V( a £p) = O (x 2 ). 
Given the linear terms in £ a , this is as good as it gets: no choice of third order terms can reduce the order of the 
Killing equation further. 

At first order, imposing the Killing identity on T requires 

Dv P aa = -(pva) + (apcr) + (aap) + d p T Vaa + Ry paa . (38) 

Note that this does not fully determine the third order terms in £ Q , but only those that are at least linear in V. To 
write the D-coefficients more explicitly, we note that at the origin of any locally inertial coordinate system (and for 
the NNCs (HI]) in particular) it holds that 

daTpjs = ~9 a dp g 7 s + d a d {l g s) p. (39) 

In addition, we may express Riemann components at the origin of a locally inertial coordinate system in terms of 
second derivatives of the metric. Substituting these expressions into (|38[) results in 

Dvpaa = -{paa) + (apa) + (aap) + d p d a g aV + -d a dv g pa - d v d {p g a ) a , (40) 
with the brackets given by 

(pea) = y~dvd v g prT + d v d ip g„ )v ^J 5%+ (^d v du g pcr - d v d (p g a)u ^J 6%. (41) 

Clearly, the bracket (paa) is symmetric in the index pair pa. It is then manifest that the RHS of (|40[) is also symmetric 
in pa, as is necessary for consistency with the symmetry of Dy paa- We remind the reader that the second derivatives 
of the metric are to be evaluated at p. These derivatives can be read off directly from the metric components listed 
in (|22]) - ((27]) and they are always proportional to a component of the Riemann tensor. Accordingly, all the coefficients 
Dvpaa are required to be a certain (sum of) Riemann component (s) evaluated at p. 

Similarly, at second order the Killing identity on Y requires the coefficients Eyvpaa to be given in terms of partial 
derivatives of the Riemann components. At the next order, coefficients Fyy y poa are required to be a sum of second 
derivatives of Riemann components and products of Riemann components, and so on. Hence we see that the Killing 
identity may be satisfied exactly on the null geodesic Y through p. This also has the consequence that £ a satisfies the 
Killing equation exactly on Y . Furthermore, as shown in [f|, £ a is null on Y and has the same proportionality to an 
affine tangent as has a Killing vector generating a Killing horizon, namely £ a |r = Vk a . 
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A. Placing the fixed point to the past of the origin 

In the treatment of [f|, the origin of the coordinate system was placed on a local causal horizon's terminal surface, 
such as the surface E we have considered here. However, it turns out to be more natural thcrmodynamically if the 
fixed point of the approximate Killing vector is placed slightly to the past of the terminal surface (see [j| for further 
details). We would therefore like to specify an approximate Killing vector that still satisfies the Killing identity on T, 
but with the fixed point now chosen to be a point on T with NNCs (0, Vo, 0, 0, ...), slightly to the past of p. Instead 
of (|34|) we expand £ a in powers of x a = x a — VoSy, as 

£ Q = U5l - V5% + \c paa x(>x° + ^D p ^ a xPx°x» + O (x 4 ) . (42) 

For simplicity we will merely demand the Killing identity on T to hold at the linear order. Then a straightforward 
calculation reveals that the .D-cocfficicnts remain unchanged, but the C-coefficients change to 

C paa = (-{paa) + (apcr) + (apa)^j V . (43) 

This does not mean however that the approximate Killing vector picks up a quadratic term, since Vq is regarded as a 
first order quantity. 
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Appendix A: A generalisation of Fermi- Walker transport 

As is well known, the inner product of a pair of vectors is preserved under parallel transport along a curve. If 
the curve is not null, any vector field defined along the curve may be uniquely decomposed into a part tangential to 
the curve and a part orthogonal to it. Since an affinc tangent to a geodesic is parallel transported into itself, this 
implies that along a non-null geodesic the tangential and orthogonal parts of a parallel transported vector are parallel 
transported separately. 

Fermi- Walker (FW) transport generalises the notion of parallel transport such that the above features survive 
along non-geodesic curves. That is, FW transport reduces to parallel transport along a geodesic, the inner product 
of a pair of FW transported vectors is preserved and the tangential and orthogonal parts of a vector that is FW 
transported along a non-null curve are FW transported separately. 

The present aim is to define a notion of transport across a multidimensional surface that again shares these features 
with FW transport. This will be achieved by generating the surface by a family of curves that emanate from a point 
and by specifying how a vector is transported along each curve. Before proceeding, it will be instructive to give a 
definition of FW transport along non-null curves 6 , as this will allow for a natural extension to the multidimensional 
case. 

Let T be the projector 7 onto a non-null curve and N = 1 — T the projector onto its orthogonal complement. A 
vector w defined along the curve may be decomposed as 

w = Tw + Nw. (Al) 

If w and x are vectors which are FW transported along the curve, their tangential and orthogonal parts should be FW 
transported separately. Therefore, the inner products of their tangential parts and of their orthogonal parts should 
both be preserved separately. For the tangential parts this reads 

= [(Tw)jTxn 

= l(Twy} a (TxT + (Tw) a {(Tx)T, (A2) 



On null geodesies, FW transport reduces to parallel transport. Non-geodesic null curves will not be considered. 
In this appendix, the abstract indices denoting a tensor will often be omitted when there is no cause for confusion. 
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where covariant differentiation along the curve was denoted by a dot. Allowing for arbitrary vectors at an initial point 
on the curve, this condition can only be satisfied provided that the covariant derivatives of the tangential parts are 
orthogonal vectors. Therefore, we impose 

(Tw)' = N (Tw)' , (A3) 

and similarly for x. The latter condition is satisfied if and only if the norm of Tw remains constant. Another statement 
equivalent to (|A3[) is that Tw is parallel transported along the curve, according to the derivative operator compatible 
with the one dimensional metric induced on the curve. Similar to the condition imposed on a tangential part Tw, we 
ensure that the inner product of orthogonal parts Nw and Nx is preserved under FW transport, by imposing that 
the covariant derivative of an orthogonal part Nw is tangential. This may be written as 

(Nw)' = T(Nw)- . (A4) 

The condition (|A3j) together with the (n — 1) conditions (|A4[) can be summarized by a single differential vector 
equation. For this purpose, we parametrize the curve such that the tangent e has a constant norm, given by e a e a = E. 
A standard choice would be to parametrize by proper length or time, for which we have E = ±1. In abstract index 
notation, the covariant derivative along the curve is defined by "." = e a V a and the tangential projector is given 
by T\ = E~ 1 e a eb- Thus the tangential part of w is written as Tw = (E~ 1 e b Wb) e. Using these expressions and 
substituting (|A3[) into (|A4[) . we find that a vector can be defined to be FW transported if it is a solution to the 
equation 

w = E- 1 [(e b w b ) e - (e b w b ) e] . (A5) 

We will now define a vector field on an m-dimensional (non-null) surface S, starting from a vector at a given point 
on the surface, and in such a way that the definition reduces to FW transport for the case m = 1. An m-dimensional 
projector T and an (n — m)-dimensional projector N are introduced, projecting onto the surface and its orthogonal 
complement respectively. The metric T a b induced on the surface is obtained from the projection operator T b a by 
lowering one index by the full spacetime metric. Furthermore, the surface is regarded as generated by a family of 
curves that emanate from a given base-point p. On a sufficiently small surface, every point (other than p) will lie on 
a unique curve in this family. We then focus on one particular curve, again denoting covariant differentiation along 
the curve by a dot. 

Starting with a vector w(p) at p, the transport of w(p) along the curve is defined by the m conditions 

(Tw)' = N(Tw)- , (A6) 

together with the (n — m) conditions 

(Nw)' = T(Nw)' ) ( A7 ) 

in close analogy to the conditions (|A3|) and (|A4[) for FW transport. That is to say, w is defined along the curve by 
requiring the covariant derivative of its tangential part to be orthogonal and the covariant derivative of its orthogonal 
part to be tangential. Condition (|A3[) may again be interpreted as stating that the tangential part of w is parallel 
transported along the curve, only now according to the derivative operator compatible with the induced metric T ab 
on the surface. By construction, inner products of transported vectors are preserved and tangential and orthogonal 
parts of a transported vector are transported separately. In particular, if w(p) is orthogonal to S then the transport 
of w(p) will be orthogonal to £ everywhere and will have a fixed norm 8 . 

Repeating this for all curves through p in the family that was chosen to generate S then obtains a vector field 
w(q) on S. Clearly, the definition of the transported vector field depends on the family of curves that generate E. 
For definiteness, we choose the family of surface geodesies through p associated to the induced metric T ab . Since a 
different choice of base-point on £ leads to a different family of surface geodesies, our definition of surface transport 
is subject to the choice of base-point p. 

As in the case of FW transport, the defining conditions (|A6|) and (|A7[) may be combined into a single equation. 
To this end, we introduce a tangential basis {e^p)}, A = 1, . . . , m at p and construct tangential basis vectors eA(q) 



In this case, it is interesting to note that the transported vector may be regarded as obtained from a continuous application of (full 
metric) parallel transport followed by orthogonal projection. This is easily seen by writing out components of the transported vector an 
infinitesimal distance away from p and comparing with the components of the orthogonal projection of the parallel transport. 
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everywhere on S as the solutions to (|A6|) with initial condition e_a(p). In other words, eA^q) is the parallel transport 
of e^G?) along the surface geodesies, according to the derivative operator compatible with the induced metric on the 
surface. As such, &a is orthogonal to E and the inner products g{eA,e.B) — Tab are constant along E. An orthonormal 
frame at p would provide a natural initial condition to give rise to such a tangential basis. Next, writing Tw = w A e^t, 
we note that the tangential components can be expressed as w A = T AB e b B Wb, where T AB is the inverse of the matrix 
Tab- Upon substitution of (|A6|) into (|A7[) . a straightforward manipulation then obtains 



For the case m = 1, this equation manifestly reduces to the defining equation (|A5|) of FW transport. Furthermore, 
(|A8[) implies the n conditions (|A6[) and (|A7[) and uniquely determines a vector w(q) along E, once w(p) is specified. 
In particular, it may be checked that the vector fields k a (q) and l a (q) defined in © are the unique solutions to (|A8[) 
for the given initial vectors k a (p) and l a (p). 
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